Large increase of the Curie temperature by orbital ordering control 
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Using first principle calculations we showed that the Curie temperature of manganites thin films 
can be increased by far more than an order of magnitude by applying appropriate strains. Our main 
breakthrough is that the control of the orbital ordering responsible for the spectacular Tc increase 
cannot be imposed by the substrate only. Indeed, the strains, first applied by the substrate, need 
to be maintained over the growth direction by the alternation of the manganite layers with another 
appropriate material. Following these theoretical findings, we synthesized such super-lattices and 
verified our theoretical predictions. 



Ferromagnetic super-lattices with very high Tc present 
obvious interests for electronic applications. This is the 
first step required to design and control magnetism in 
more elaborate structures. Perovskite manganese ox- 
ides of generic formula Ai_xBxMn03 (A being a triva- 
lent cation and B a divalent one) are metallic and fer- 
romagnetic in the approximative range of doping values 
0.17 < X < 0.5 [5]. These properties are directly related 
to the 3d shell of the manganese atoms and can be easily 
explained by the double-exchange mechanism. The Curie 
temperature is thus directly related to the delocalization 
within the Cg orbitals of the Mn^+ / Mn'^+ ions [2], that 
is to the manganese ions orbital order. Finding a way 
to design such materials with the desired 3d orbital or- 
der would thus open the way to the design of artificial 
materials with the desired magnetic properties, and thus 
to novel magneto-electronic applications, including spin- 
valve devices or non- volatile magnetic memory working 
far above room temperature. 

The manganese Cg shell is doubly degenerated (in a 
regular octahedron environment) however only partly oc- 
cupied. The Mn^+ ions are thus Jahn- Teller active and 
small distortions of the oxygen octahedra can stabilize 
one of the Cg orbitals with respect to the other. If one 
can impose a specific geometry around the Mn ions (for 
instance using synthesis in thin films), it should thus be 
possible to choose between the d^z'^-r'^ and the dx2_y2 
orbitals which one will be occupied and which one will 
be empty. If the d^z'^-r'^ is occupied, the double exchange 
will take place essentially along the c direction, that is 
between (a, 6) planes (see figure [1]). On the contrary if 
the dx2_y2 is stabilized, the intra-layer double exchange 
will be very strong and the inter-layer weak. In an es- 
sentially 2-dimensional (2-D) system such as thin films, 
only the in-plane interactions are important and the con- 
trol of the orbital ordering should thus allow the control 
of magnetic exchange and the associated Curie temper- 
ature, Tc. According to the above argument, the Curie 
temperature of 2-D thin films will be maximized if one 



(i) stabilizes the dx2_y2 orbital over the d^z'^-r'^ and (ii) 
maximizes the in-plane effective exchange integral, be- 
tween the dx2_y2 orbitals, by preventing the tilt of the 
octahedra observed in the bulk 0. Indeed, in the double 
exchange mechanism, Tc scales as 
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(J: effective exchange between adjacent Mn ions, tdp'. 
{Mn)d-{0)2p transfer integral and Sdp- overlap) and 
the octahedra tilt lowers the Sdp overlap between the 



Mn d^ 



and the O 2px bridging orbitals. 



(oa) 



Weak in plane delocalisation 
Weak Curie temperature 



(c<a) 




d^2_y2 



X_L_L 




* Deceased february 2009 



Large in plane delocalisation 
Large Curie temperature 



FIG. 1: (Color online) Schematic representation of orbital or- 
dering induced by cell distortions in a two dimensional system, 
(c > a) The elongation in the c direction favors the d^z'^-r^ 
orbital occupation and weak in-plane ferromagnetic coupling, 
the delocalisation being essentially in the c direction, (c < a) 
The contraction of the c parameter favors the d^2_y2 orbital 
occupation, strong in-plane ferromagnetic coupling and large 
Curie temperature. 



Ai-xBxMnOa thin films deposited on a cubic substrate 
presenting a = h parameters slightly larger than the bulk 
manganite value should thus insure an increase of the 
manganite in-plane parameters and, in order to approx- 
imately conserve the cell volume, a decrease of the c pa- 
rameter. With c/a < 1 the dx2_y2 orbital occupation 
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and large Tc should be favored, as desired. 

Several authors synthesized such thin films as for in- 
stance La2/3Sri/3Mn03 (LSMO) thin films deposited on 
the SrTiOa cubic substrate They found that for 

a small number of LSMO unit cells in the c direction, 
the magnetization as well as the Curie Temperature are 
strongly decreased compared to the 370K bulk value 
From twelve deposited unit cells however the bulk be- 
havior is approximately recovered with just a slight Tc 
reduction. It is clear that the expected behavior and 
thus orbital ordering is not realized in these films. In 
fact the quick recovering of the bulk magnetic behavior 
tells us that the strains imposed by the substrate are es- 
sentially very quickly relaxed as a function of the film 
thickness. This is most probably the case for the octa- 
hedra tilt present in the LSMO bulk. One should thus 
find a way to maintain the desired strains (in-place cell 
parameter elongation and quadratic symmetry) over the 
whole film. For this purpose one can think to alternate 
LSMO layers with layers of an adequate cubic material as 
in [(La2/3Sri/3Mn03)n(BaTi03)p]^ super-lattices. We 
thus explored this possibility using first principle calcu- 
lations. 

We studied, using periodic density functional the- 
ory the (La2/3Sri/3Mn03)3(BaTi03)3 super-lattice 
deposited on a lattice matched SrTi03 (OOl)-oriented 
substrate (see figure [2]). Since the epitaxial films nor- 
mally follow the structure of the perovskite substrate, 
we imposed to our system the substrate in-plane lattice 
constants, a = b = 3.9056 A. The c parameter was opti- 
mized as well as atomic positions. 

The super-lattice is built from metallic LSMO layers 
and insulating BaTi03 (BTO) layers. It is thus of cru- 
cial importance for the reliability of the results to be able 
to accurately position the LSMO Fermi level referring to 
the BTO gap. For this purpose we used hybrid function- 
als known to correctly reproduce insulating _gaps, that 
is B3LYP [13] and the very recent BIWC Q specifi- 
cally designed to obtain the correct properties on fer- 
roelectric materials such as BaTi03. We thus used the 
CRYSTAL package ^ with the basis sets and effective 
core pseudo-potentials (ECP) of reference 15. Another 
important technical point is the treatment of the La/Sr 
disorder. Indeed, in periodic calculations it is not possi- 
ble to properly mimic the ionic disorder. We thus made 
a whole set of calculations associated with the two ex- 
treme cases, that is with the Sr ions totally ordered (3 
different positions) and with average La/Sr ions. These 
average ions were modeled either using ordered La/Sr 
ECPs but with averaged effective nuclear charges (3 pos- 
sibilities) or by using only La ECPs with average effective 
charges. All 7 different approximations of the ionic dis- 
order were computed with the two functionals. All of 
them yielded equivalent results and physical properties, 
even if one finds slight numerical differences. One set 
of numerical values (corresponding to B3LYP functional. 



ordered ECPs with the strontium in the middle layer and 
averaged charges) are presented as a matter of example. 
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FIG. 2: (Color Online) Schematic representation of the super- 
lattice and calculated inter-planar distances, c', and associ- 
ated sub-cell volumes, V for the LSMO layer (A). 

Considering the fact that the LSMO layers are ex- 
tended in the (a, 6) directions compared to the bulk 
{a = b = c = 3.87 A), our results (figure [2]) are at first 
quite surprising. Indeed, the two LSMO layers at the in- 
terfaces with the BaTiOs are not only extended in the a 
and b directions, but also in the c direction. It results a 
strong extension of the associated sub-cell volume with 
respectively 66.57 A^ — for the interface with a Ba-0 
layer — and 61.44 A^ — for the interface with a Ti-0 
layer — to be compared with the bulk pseudo-cubic unit 
cell volume of 58.3 A^ 0]. On the contrary, the central 
LSMO layer is strongly contracted in the c direction re- 
sulting in a reduction of the associated sub-cell volume 
(56.87 A^). One sees immediately that for the interface 
layers the c'/a > 1 anisotropy strongly favors the stabi- 
lization of the 6/3^2 orbital. This stabilization is indeed 
associated in our calculations with a larger occupation 
of the d^z'^-r'^ orbital compared to the dx2-y2 one. The 
anisotropy is reversed for the central layer {c' /a < 1) and 
thus the dx2_y2 orbital is stabilized and presents a larger 
occupation number. The magnetism is essentially two- 
dimensional, supported by the Mn ions (with 3.467//^, 
3.541/i^, 3.886/iB)- The BaTiOa layers remain essen- 
tially non magnetic. We can thus predict two magnetic 
behaviors in these super-lattices. The first one is related 
to the interfacial layers and is associated with a weak 
in-plane double-exchange and thus a weak Curie's tem- 
perature (Tci). The second one is related to the central 
layers and is associated with a strong in-plane double- 
exchange and thus a strong Curie's temperature (Tc2). 
One can thus predict three different phases, the first one 
at very low temperature with all three LSMO layers par- 
ticipating to the magnetization, the second one at larger 
temperature with only the central layer participating to 
the magnetism and finally, at high temperature, the para- 
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magnetic phase. 

In an attempt to estimate the order of magnitude 
of Tci/Tc2, one can estimate the double-exchange mag- 
netic coupling. Indeed, even if in bulk manganites, the 
double-exchange is not always sufficient to explain the 
physics and polaron formation should be considered [l^ , 
in fully- strained films such as the present one, we do not 
expect such relaxation mechanism to take place. The 
double-exchange integral is directly proportional to the 
delocalization between adjacent magnetic ions and thus 
to the square of the overlap integral between the mag- 
netic 3d orbital and the ligand 2p bridging orbital. The 
Tci/Tc2 ^ \{Mnds,2_^2\0 2p^)/{Mnd^2_y20 2p^)f ratio 
can thus be undervalued by l/(4>/3)^ = 1/48. Of course, 
this is a calculation using crude approximations, however 
it tells us that there is more than an order of magnitude 
between Td and Tc2. 

Another important result is the fact that the BaTiOa 
triple-layer does not preserve its insulating character but 
exhibit a weak, however non null, density of states at the 
Fermi level in the dominant LSMO spin orientation. This 
electronic delocalization through the BaTiOs triple-layer 
insures the ferromagnetic coupling between the different 
LSMO layers by what we could call a long-range double- 
exchange mechanism. 

As a conclusion, ffist principles calculations predict for 
the super-lattice 

• two different ferromagnetic phases, associated with 
two Curie temperatures differing by nearly two or- 
ders of magnitude, 

• an optimized high Curie temperature expected to 
be larger than the bulk one (no octahedra tilt). 

We thus synthesized (La2/3Sri/3Mn03)3(BaTi03)3 
super-lattices using laser-MBE deposition [17] on a [001]- 
oriented SrTi03 substrate. We also synthetised similar 
super-lattices using the isoelectronic calcium manganite: 
La2/3Cai/3Mn03 (LCMO). Indeed, one can expect in the 
LCMO super-lattice similar orbital ordering effects, for 
similar reasons, as predicted in the strontium one. 

The pulsed laser deposition was performed in low 
pressure (5.10~^mbar) to enable control of the growth 
using the Reflexion High Energy Electron Diffraction 
(RHEED). While the growth of the structures was car- 
ried out in low ozone content (0.1% volume), to promote 
a maximum oxidation of the super-lattice without fa- 
vor izing a strong diffusion between the different layers, 
the structures were rapidly cooled in high ozone con- 
centration (7% volume). The deposition of each mate- 
rial, namely BTO and LSMO or LCMO, was measured 
and calibrated by growing them independently under the 
same conditions (620°C, pressure 5% x 10~^mbar, ozone 
concentration 0.1% volume) using an in situ RHEED 
control for a precise control of the growth of each layer. 
Figure [3] shows two RHEED patterns, taken along the 



[110] azimuth of the substrate. These two pictures evi- 
dence a streaky-like diffiaction pattern in agreement with 
a 2D growth of the super-lattice. The comparison be- 




FIG. 3: RHEED patterns for the BTO/LSMO super- lattice 
along the [110] azimuth of the SrTiOa substrate. A: 620° C, 
5 X 10"^mbar (O2 + O3). B: 120°C, 5 x 10"^mbar. 

tween the two patterns taken at the end of the deposi- 
tion (A) and at the end of the cooling (B) evidences two 
features. First, there is no drastic modification of the 
pattern after the cooling process. Second the in-plane 
structure of the super-lattice exhibits the same lattice 
parameters as the SrTi03 substrate and thus can be de- 
scribed in a perovskite-like sub-cell. The B pattern also 
presents Kikuchi lines indicating the high quality of the 
deposited film. These lines are not observed in pattern 
A, due to the important scattering of the electrons at the 
deposition pressure. 

To achieve a better structural characterization of 
the super-lattices an X-ray diffraction study was car- 
ried out using a 4-circle diffractometer. Two dif- 
ferent set of reflections were studied: the out of 
plane reflections and the so-called asymmetric reflec- 
tions. From the first set of reflections, the value of the 
out-of-plane lattice parameter could be calculated for 
the [(La2/3Sri/3Mn03)3(BaTi03)3]25 super-lattice. We 
measured an out-of-plane lattice parameter of 24.281 A, 
corresponding almost (1%) to the value determined from 
the theoretical calculations (24.519 A). This value is 
larger than the expected value, calculated from the av- 
erage out-of-plane parameters of films of the two start- 
ing materials. A Wiliamson-Hall analysis carried out on 
the [001] reflections shows that the coherence length of 
the diffraction along that direction corresponds to 600 A, 
that is the thickness of the deposited film. The same 
analysis carried out on asymmetrical reflections ([Oil], 
[022], ...) indicates that the in-plane lattice parameter is 
3.9048 A, in agreement with the substrate lattice parame- 
ter, and that the coherence length of the diffraction along 
an in-plane direction is 1000 A, in agreement with the size 
of the terraces observed on the SrTi03 substrates. 

The magnetic measurements (figure H]) , carried out in 
the temperature range 4-1000 K, evidence two ferromag- 
netic phases, as predicted. The first one at low temper- 
ature and the second for Td < T < Tc2. The low Curie 
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temperature, expected to correspond to the two interfa- 
cial manganite layers, is about 25 K for the LSMO 

super-lattice and 50K for the LCMO one. The 

high Curie temperature, expected to correspond to the 
central manganite layer, is about 1000 K for the LCMO 
system and about 650 K for the LSMO one, that is far 
above both Td and the bulk Curie temperatures, as pre- 
dicted by theoretical calculations. Indeed^he Curie tem- 
perature of LSMO bulk or thin films H is 370 K and 
in the range 250-260 K for LCMO bulk [1, [l^ and thin 
films The two inserts of figure [4] show that magnetic 
loops can be registered both at 10 K and 300 K, insuring 
true ferromagnetic phases. 
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FIG. 4: (Color Online) Magnetization (500 Oerst- 
eds) for the [(La2/3Sri/3Mn03)3(BaTi03)3]25 (red squares) 
and [(La2/3Cai/3MnO3)4(BaTiO3)4]20 (green circles) super- 
lattices. Hysteresis loops of the calcium super-lattice, 
recorded at lOK and 300K, are presented in the inserts. The 
arrows mark the bulk Curie temperatures. The surface of 
both samples is 16mm^. 

In summary we showed that it is possible to increase 
the Curie temperature of a manganite ferromagnetic thin 
film far above the bulk Curie temperature, by the control 
of orbital ordering. The key idea to reach this result is 
that the strains applied by the substrate are not suffi- 
cient by themselves and that they must be maintained 
using the regular intercalation of adequate layers. Our 
findings open new opportunities for the design and con- 
trol of magnetism in artificial structures and pave the 
way to novel magneto-electronic applications, including 
non-volatile magnetic memory working far above room 
temperature. 
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